Straining the identity of Majorana fermions 
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We propose an experimental setup of an interferometer for the observation of neutral Majorana 
fermions on topological insulator — superconductor — ferromagnet junctions. We show that the 
extended lattice defects naturally present in materials, dislocations, induce spin currents on the 
edges while keeping the bulk time-reversal symmetry intact. We propose a simple two terminal 
conductance measurement in an interferometer formed by two edge point contacts, which reveals 
the nature of Majorana states through the effect of dislocations. The zero temperature magneto- 
conductance changes from even oscillations with period </>o/2 (0o is the flux quantum hc/e) to odd 
oscillations with period <f>o, when non-trivial dislocations are present and the Majorana states are 
sufficiently strongly coupled. Additionally, the conductance acquires a notable asymmetry as a 
function of the incident electron energy, due to the topological influence of the dislocations, while 
resonances appear at the coupling energy of Majorana states. 
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There is a strong interest to realize, observe and ma- 
nipulate Majorana fermions, because of the non-Abclian 
statistics they possess |l|, being the basis for topologi- 
cal quantum computation [2j. Majorana fermions have 
been argued to be present in the v = 5/2 fractional 
quantum Hall state [l|, |3(] , in the p— wave superconductor 
Sr2Ru04 14 and in topological insulator-superconductor 
junctions p, Q . Topological insulators (TI) [7| have gap- 
less edge (2DTI) or surface (3DTI) states that are heli- 
cal and topologically protected in the absence of time- 
reversal symmetry (TRS) breaking fields. Breaking TRS 
by depositing an insulating magnetic (M) material can 
open an energy gap leading to a novel surface quantum 
Hall effect with a xy = ±e 2 /2h [8|. Further, deposition of 
a superconductor (S) on the edge or surface leads to Ma- 
jorana bound states (MBS) at an S-TI-M interface where 
the gap changes sign [5j . These Majorana fermions have 
proven to be very elusive since they are neutral, and there 
are a few proposals for their observation, ranging from 
rather indirect tunneling experiments |9j, to interference 
experiments [1CH12J . In this paper, we make a step fur- 
ther and propose a novel interferometer which preserves 
TRS and can be used not only for identifying the Majo- 
ranas, but also for the dual purpose of understanding the 
fundamental properties of topological lattice defects. 

We propose a standard Aharonov-Bohm (AB) inter- 
ferometer (Fig. [TJd), where the presence of dislocations 
within the interferometer area causes a topological phase 
shift on the edge states due to the translational effect of 
the dislocation Burgers vector on the edge wavefunction. 
This AB effect I13J is analogous to the effect of pierced 
magnetic flux [lOj, except that it preserves TRS. The 
magnetic flux induces electrical current flow, the persis- 
tent current, in the ground state. Analogously, disloca- 
tions induce the TRS invariant counterpart, dissipation- 
less spin currents. Spin currents are typically hard to ob- 
serve, but appear useful for MBS detection. Dislocations 



in 3DTI were also found to host interesting states [14[. 

The STIM interface locally breaks TRS and particle- 
hole symmetry (PHS) |5j, |6|, so that clear experimen- 
tal signatures in the two-terminal AB interferometer are 
expected, e.g. asymmetry of the magnetoconductance 
(G(4>) y^ G{— </>), where <p is the threaded magnetic flux), 
being typically absent due to TRS. We find that mag- 
netoconductance remains even in the presence of MBS, 
due to the topological helicity symmetry (exchange of the 
left/right moving up/down spin, for the left/right moving 
down/up spin edge modes). But when dislocations are 
present (this is controlled by straining the bulk of the TI, 
Fig UtO , a spin current is introduced in the interferome- 
ter which is sensitive to the helicity flip and therefore can 
detect the signatures of MBS. Most strikingly, the oscil- 
lations 5G(4>) switch from even with period cj>o/2 (0o is 
the flux quantum hc/e) to odd oscillations with period 4>o 
when dislocations enter the device and the MBS are cou- 
pled (Fig. [3]), while oscillations vanish in absence of MBS 
at the STIM. We predict that the conductance satisfies 



G(<l>,E,<l> d )=G(-<f>,-E,-<l> d ) 



(1) 



(E incident electron energy, and 4>d the dislocation scat- 
tering topological phase), which allows the use of the 
topological effect of dislocations (<f>d) as a new control 
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parameter (absent in all existing proposals 
bring out the signatures of the MBS. 

In the example of 2DTI realized in Hg(Cd)Tc quantum 
wells [71 Il5|. dislocations seem neatly controllable. Ac- 
cording to detailed structure studies as in Ref. [16| , at low 
temperatures and at the yield stress of ~ 10 — lOOMPa, 
there are 10~ 10 m~ 2 dislocations, giving < 1 mobile dislo- 
cations piercing a 5/xm x 5/mi sample. After yielding, no 
additional stress is needed in Hg(Cd)Te to move defects, 
so they move freely and independently Upon reaching a 
high 10% strain, with a total Burgers vector of 2 * 10 _9 m, 
one expects that ~ 150 dislocations have passed (glided) 
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FIG. 1. Experimental setup for observing neutral Majorana 
fermion bound states (MBS), a) Three-point bending (by 
force F), moving an example dislocation line (Burgers vector 
b) which pierces the substrate (green) and the 2D Topological 
Insulator (TI, grey) at the circle. At one TI edge there is a 
pair of superconductor (SC, yellow strip) - ferromagnet (FM, 
orange strip) - TI junctions. Contacts are marked violet, b) 
Edge modes of the 2DTI (grey area) traverse the interfer- 
ometer, where the MBS (red dots) are present in the upper 
arm. The dislocation induces a translation on the half-plane 
of missing atoms (dotted line), causing a topological phase 
shift exp (iK ■ b), where K is the three-dimensional embed- 
ding of the edge Fermi momentum K. Inset: Schematic 2DTI 
band structure showing K near n. 



through such a sample. With dislocations being the most 
natural and abundant topological defects in crystals, we 
further expect the rightful use of dislocation induced spin 
currents as novel TRS probes in the future. 

Our interferometer is made of a 2DTI shaped by two 
point contacts (Fig. [Tb), and we model it using the 
Landaucr-Buttiker scattering matrix formalism valid at 
low temperatures in the regime of coherent transport [17| . 
Edge segments comprising the interferometer support 
one electron and one hole chirai mode. The Bogoliubov - 
de Gennes Hamiltonian describing each edge segment is, 

T 3 (v F p a 3 + A d + T 3 a 3 eA/hc - E F )^ = EV, (2) 

where p = —ihd/dx, Ep the Fermi energy, vf the 
Fermi velocity, and A the magnetic vector potential, and 
the x-axis is along the given edge segment. The four- 
component spinor is ^ = (\Pei-; ^ei, &hft ^h\) T while the 
r matrices mix the electron and hole parts of the wave- 
function, and a the spin components. 

The effect of dislocations is contained in the po- 
tential Ad of Eq. ([2]). It encodes for the AB effect 
exp(^- § Ad dx) — exp (i2it(j)d), with pseudo-flux (f>d 
stemming from the topological effect of the dislocation 
on the wavefunctions on the edge. It is well known that 



this effect is described by a translation by the Burg- 
ers vector b on traversal of electron around the dislo- 
cation core line threading the TI inside the ring-shaped 
area of the interferometer 18( . The translation operator 



exp (iK ■ b) = exp (i2n(j)d) is determined by the three- 
dimensional Burgers vector b of the dislocation line that 
can be of any type (edge, screw or mixed). The vec- 
tor K is the three-dimensional embedding of the edge 
wavefunction wavevector K [7[ (see Fig. QJi). The dis- 
location effects discussed in this paper depend on one- 
dimensional momentum K on the edge being non-zero; 
such 2DTI variety is not yet observed, but it could exist 
in Hg(Cd)Te [19J], or Heusler alloys [20(. Dislocations pre- 
serve time reversal and particle hole (PHS) symmetries, 
represented by T = ia^C and 5 = T2CT2C, respectively, 
with C the complex conjugation, and they are distinct 
from ordinary disorder due to their intrinsic gauge sym- 
metry. Gcncrically, as in HgTc wells, edge segments do 
not exhibit PHS, so we checked that our results are ro- 
bust to breaking PHS by assigning different velocities to 
edge states below and above E = (e.g. the edge energy 



spectrum observed in 3DTI |21|). 

Scattering formalism- The Hamiltonian of Eq. (J2j) de- 
termines the energy dependent wavevector of the left 
(spin down on upper edge), and right (spin up on upper 
edge), moving electron, as well as their time reversed hole 
pairs. The point contacts, the two halves of the upper 
ring arm, the coupled MBS between the two upper arm 
halves, and the lower ring arm are all described by single 
scattering points with corresponding matrices (S S catt)- 
Each matrix S sca ,tt connects the amplitudes (Oj, and 
On) of the modes outgoing to the left/right (L/R) side, 
to the amplitudes (II and In) of the incoming modes, 
with respect to that particular scatterer. Using T = 
(0 Ll O R ) T = (<&,<&, <#,<&) T and I T = (I L ,I R ) T 
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r t , \ J , representing re- 
flection (r, r') and transmission (t, t'), where each block 

has electron/hole components, e.g. t = ( *^ l h l ) . When 
TRS is obeyed, backscattering is forbidden in the single- 
particle formalism, r = r' = [22J . We compute the total 
scattering matrix S for the four leads (labeled from 1 to 
4, Fig. [Tb), determining the conductance of the device. 

Particle conservation is enforced by S^S = 1. For 
scattering matrices connecting two edge segments, TRS 
demands S , sca tt(0) = -a 3 S sca tt(~4>) T tt 3 , and PHS 
S sca tt(E) = PiSecott(—E)*Pi, where a and /3 arc Pauli 
matrices acting on the L/R and e/h indices of S r SC att, re- 
spectively. For scattering involving all four edges (like 
in S) one should only replace a by a <£> a' , where a 1 
matrices exchange the two leads on the same side (i.e. 
1 and 4, or 2 and 3). The scattering caused by the 
coupling to, and propagation through the two MBS in 
the upper arm is given by the scattering matrix Smbs 
found in Ref. [6(. It is determined by two energy scales, 
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FIG. 2. The energy (voltage) dependent conductance, as func- 
tion of dislocation, with Em = 1. (a) Typical curves: no dis- 
location (dashed thin line); dislocation phase (j>d = 0.1 (thick 
black); and <f>d = —0.1 (thick light grey (blue)). Note the res- 
onances at izEm = ±1, and the dislocation induced asymme- 
try, (b) The asymmetry of the G(E) curves (full black line), 
calculated as Max E [G(E) - G(-E)] on the interval E G [0, 3] 
(it reaches 2 for purely odd G ~ sin(E)). The dashed grey 
(blue) line shows the amplitude of G(E) oscillations around 
the mean. The curves are robust to changes in Em- 



the coupling between the two MBS Em, and the cou- 
pling of edge states to the MBS I\ Length is mea- 
sured in units of the ring circumference L, <f> in units 
of the flux quantum 4>o = hc/e, and energy in units of 
Hvf/L. Wc consider the scattering mechanisms as fol- 
lows: (a) Propagation in the lower arm S\ ow is determined 
by nonzero elements t^ w = cxp [ild(E + 2w(j)d — 2-7T0)], 



hi, 

low 



exp [ild(E — 2-K(j)d + 2ir </>)], where l u is the length 
of the lower arm; (b) In the upper arm segments S up = 
Siow> with Id replaced by l u \, and l u n in the two segments, 
respectively; (c) Without loss of generality we take the 
point contact scattering matrix Spc to be real and satis- 
fying the TRS and PHS symmetries (edge segments are 
ordered as (1,4) on the left and (2, 3) on right, cf. Fig.QJ: 



Spc = 



/ a b b'\ 

a -b' b 

b b' -a 

\-V b -a / 



(3) 



with j3o — 1, and a 2 + b 2 + b' 2 = 1. Parameter a de- 
scribes the coupling of the ring-shaped middle of the in- 
terferometer to the leads (a = corresponds to an iso- 
lated ring with G = 0). The ratio e = b/b' measures the 
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FIG. 3. Magnetoconductance G(cj>) at zero energy (voltage) 
as function of dislocations and Majorana couplings. In ab- 
sence of MBS, G(4>) vanishes (dotted black line). The dash- 
dotted (red), dashed (blue) and full (black) curves correspond 
to three regimes of MBS coupling Em — 0,0.1,0.3, respec- 
tively, with respect to MBS-edge coupling r = 0.1 (units 
Hvf/L). The case of absence (cpd = 0) or presence (<f>d = 0.3) 
of dislocations is distinguished by thin and thick lines, respec- 
tively, for each Em value. Without dislocations, the result 
5G ~ cos (20) is robust to interferometer parameter changes. 
The presence of dislocations affects only the cases of coupled 
MBS, Em 7^ 0, by suppressing the 0o/2 harmonic in G(4>) in 
favor of the <f>o, which is always odd, i.e. sin (</>). 



asymmetry of current injected into the lower and upper 
ring arms (e = corresponds to all particles from lead 
1 being injected into the lower arm, and all from lead 
4 into the upper). Following Refs. [23j, in the present 
single particle scattering we attain the conductance of 
the charge-conductor/spin-insulator (CI) state by choos- 
ing a = l/v3, e = 1) being in the regime of Luttingcr 
liquid coupling g c > 2. In the realistic case of interme- 
diate < a < 1, the dependence on a and e is weak, 
so we present results for CI contacts. The conductance 
is given by G = e a /hY,*^.* {\Sf? | 2 - |S#f ), where i,j 
label the leads, and holes contribute opposite charge cur- 
rent from electrons. The zero temperature conductance 
at zero voltage corresponds to taking E = 0, while at low 
temperature and voltage difference, E is given by the ex- 
ternal voltage (E = eVi). We consider Ep = (small 
Ep is negligible when K ~ 7T, see inset Fig. Q]), and fix 
lui = lu2 = ld/% = L/2, while the results are insensitive 
to the asymmetry in l u \ and l u i- The point contact pa- 
rameters a, e are set to be the same in the left and right 
contact, since results are robust to this asymmetry too. 

The symmetry expressed in Eq. ([1]) is most revealing 
since it controls the behavior of the conductance G(<p, E), 
given the changes in the net Burgers vector d. It repre- 
sents the invariance of the edge states to switching the 
spin orientation of left and right moving carriers. (This 
orientation is set by the sign of the bulk spin-orbit cou- 
pling.) For the scattering on the edge, this switch is 



represented by conjugation C, i.e. the combined time- 
reversal and spin-flip operation. In this case, it follows 
that S sca .tt(E) = Sscatt (-E)* and non-trivially for the 



case of SmbSi this property holds because H^ 
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-H 
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The two-level Hamiltonian Hm fundamentally obeys the 
relation because the Majorana fields are real (7^ = j a ). 
We expect the spin-flip symmetry to be robust in absence 
of Zeeman type coupling to out-of-plane magnetic fields. 

We first consider the effect of dislocations on a triv- 
ial interferometer, one without a STIM interface. The 
presence of ^ ^ introduces a deviation from even- 
ness in G(E), as the symmetry G(E, <fi d ) = G(—E, —<j>d) 
suggests. The magnctoconductancc G{4>) = G(— 0) 
stays even, protected by TRS in a two terminal mea- 
surement. However, the nature of the G(0) oscillations 
switches from dominantly universal conductance fluctua- 
tions (UCF), i.e. period 0o, to a dominantly period 0o/2 
nature, when dislocation is introduced. 

Secondly, we introduce the STIM interface into the 
upper arm of the interferometer (cf. Fig. Q]). If there are 
no MBS forming, the STIM is a segment of gapped edge 
states with a TRS violation. The absence of MBS is 
modeled by setting T = (decoupling from the edges). 
In this case, the magnctoconductancc oscillations 
8G(4>) vanish. The dislocations influence the oscilla- 
tions, and G(E,<t> d ) = f(E - d ), with f{x) = b' 2 (l - 
a 4 ) [(l + a 2 )(l + e 2 )+4aecos(x/2)] / [l + a 8 + 2a 4 cos (2a;)] 
shows clearly that the asymmetry of G(E) is controlled 
by the dislocations. The effect persists in the limit where 
the central ring is decoupled from the leads (a = 0): 
the spectrum of the ring is given by the solutions 
of cos (2E) = cos (ir(j>d), and shows the symmetries 
{E n {4> d )} jL ~{E n {cp d )} and {E n {cf> d )} = -{E n (-<f> d )}. 

The general asymmetry features in G[E) due to dislo- 
cations persist when MBS are added, and new signature 
effects appear in the magnctoconductance as dislocations 
are manipulated. The G(E) shows oscillatory behav- 
ior, with resonances at ±£jv/, shown in Fig. Ela). In 
Fig. HJb), we provide a summary of the dislocation effect 
on the behavior of G(E). Introduction of non-zero dislo- 
cation phase causes a large asymmetry that persists for 
all values of Em- If the flux is present, G(E) becomes 
asymmetric at any value of <p d , and more strongly as Em 
increases (note that when MBS are absent, there is no 
dependence on 0). The last observation was made also 
for a more complicated hypothetical interferometer [11] . 

Fig. [3] presents the characteristic influence of disloca- 
tions and Majorana states on the magnctoconductance 
at zero energy (i.e. zero voltage at low temperatures). 
As announced, even though TRS is broken by the MBS 
scattering, a resulting non-even G(0) is observed only 
in the presence of dislocations. Namely, SG(<fi) has two 
prominent Fourier components, and both have definite 
parity: the UCF in the form of sin(0), and the har- 
monic cos (20). When Em = (MBS decoupled from 
each other), the UCF vanish. However, when Em 7^ 0, 



dislocations show a clear signature: in their presence, as 
Em /r increases, the harmonic is suppressed in favor of 
the UCF, and therefore simultaneously the transforma- 
tion from even to odd G(0) is observed. If E M > T, a 
small value of cf> d (e.g. 0.05) already causes a linear G(0) 
up to |0| < 1/4 (cf. Fig.©. 

In conclusion, we demonstrated the usefulness of shear 
stress manipulated dislocations in observing neutral Ma- 
jorana fermions in TIs. We found clear signatures of 
dislocation - MBS interplay in magnctoconductance os- 
cillations at zero energy, and showed the enhanced con- 
ductance symmetry of Eq. ([T]), a direct consequence of 
the symmetry of the TI and the presence of dislocations. 
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